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Abstract 

We show that the gauge- fermion interaction in multiflavour (2 + l)-dimensional quantum 
electrodynamics with a finite infrared cut-off is responsible for non-fermi liquid behaviour in the 
infrared, in the sense of leading to the existence of a non-trivial fixed point at zero momentum, 
as well as to a significant slowing down of the running of the coupling at intermediate scales as 
compared with previous analyses on the subject. Both these features constitute deviations from 
fermi-liquid theory. Our discussion is based on the leading- 1/A^ resummed solution for the 
wave-function renormalization of the Schwinger- Dyson equations . The present work completes 
and confirms the expectations of an earlier work by two of the authors (I.J.R.A. and N.E.M.) 
on the non-trivial infrared structure of the theory. 
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Recently two of us in a paper from now on to be referred as I, have suggested that a 
Schwinger-Dyson (SD) analysis for the self-energy of the massless fermion of QED3 provided 
intriguing evidence for a non-trivial infrared structure of the theory. We argued that the effec- 
tive dimensionless running coupling constant varied slowly over a wide range of intermediate 
momenta lying between a (the mass scale of Q^^-Ds, which plays the role of an ultra-violet cut- 
off 0) and the infrared cut-off e. This situation resembles that in (four-dimensional) 'walking 
technicolour' models of particle physics, and we described it in terms of a 'quasi-fixed point'. 
The inclusion of wave-function renormalization was crucial to the effect. Equally important 
was the incorporation of a finite infrared cut-off. 



This situation is very interesting from a condensed-matter point of view. There have been ar- 
guments that a variant of QED3 may serve as a continuum model for simulating the physics 
of the novel high-temperature superconductors. On the one hand, dynamical mass generation 
provides a mechanism for the spontaneous breaking of electromagnetic U{1) symmetry, thus 
leading to superconductivity. On the other , deviations from trivial infrared fixed point struc- 
ture have recently been argued to be responsible for departures from non-fermi liquid behaviour 
in the normal phase 0, |l| in which no dynamical mass is generated. Indeed the presence of 
wave-function renormalization yields marginal 0(l/A^)-corrections to the 'bulk' non-fermi liq- 
uid behaviour caused by the gauge interaction in the limit of infinite flavour number ^ 00. 
Such corrections lead to the appearance of modified critical exponents. In particular, at low 
temperatures there appear to be logarithmic scaling violations of the linear resistivity of the 
system of order 0{1/N) |]T|. 

In I we based our arguments on an approximate form of the leading- resummed SD 
equations, which was an improved version of the results of ref. 0. The approximation employed 
in 1^ was strictly valid for a regime of momenta p << a, while that in was appropriate to 
p ~ a. For the actual physical models, however, the regime of momenta over which deviations 
from fermi liquid behaviour is observed is much wider than those employed in ref. |^ or ref.|^]. 
To extend the validity of the arguments of I over a wider momentum regime in a mathematically 
rigorous way one needs to have a more precise picture of the situation by employing as exact 
expressions for the 0{1/N) SD equations as possible. In particular we would like to re-examine 
the way the effective coupling constant runs in the momentum range < p < a, in a regime of 
parameters such that there is no dynamical mass generation. That is the purpose of this note. 
As we shall see, the situation is exactly as predicted by the preliminary approximate analysis 
of I, thereby confirming the expectations of that reference. 



To this end, it seems instructive first to recapitulate the arguments of I while briefly reviewing 
the formalism. The SD equations, in the one-loop resummed limit, read 
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where the Landau gauge has been assumed, and a = is the dynamically-generated scale 
of the theory. A{p) is the wave- function renormalization, B{p) is the mass-gap function, and 
k"^) is a vertex function. The integrals in (|l|) and (0) are effectively cut-off at a, due to a 
sharp decay of the integrands above this scale [|]. Hence a may be considered as an effective 
UV cut-off. Moreover, for reasons that will become clear below we also introduce an infrared 
cut-off e. Following ref. 0, we make the vertex ansatz 

T,{p',k')=^,A{kr^j,G{k'). (4) 

The Pennington and Webb ansatz corresponds to = 1, where chiral symmetry breaking 
occurs for arbitrarily large iV [0. It is this case that was argued to be consistent with the Ward 
identities that follow from gauge invariance 0. Here we shall concentrate on the n = 1 ansatz. 
It is interesting to note from (|ip,(^) that in this case the right hand sides of the SD equations 
depend only on the ratio B{p)/A{p) = which is the physical gap function. 

For the purposes of I, and of this note, we are interested in the normal phase of the theory 
where the mass gap B{p), or equivalently can be ignored in the denominators of (|T]). In 

that case,equation (|l]) for A{p) reduces to simple quadrature if n = 1. Using the ansatz (^), 
one can then analyze the Schwinger- Dyson (SD) equations, in the various regimes of momenta, 
in terms of a running coupling obtained from substituting the solution for A{p) into equation 
(I) for the gap: 

9r{p, e) = -JT^ 
A{p, e) 

where Qq = S/n^N, N is the number of fermion flavours, and e is an infrared cutoff. The reader 
might have wondered whether the definition (|^) of a running coupling makes any sense in the 
normal phase of the model, where dynamical mass generation is absent. The resolution of this 
comes from treating the gap equation as a non-trivial equation, of which the vanishing gap 
appears as a consistent solution. This provides a definition of gji which continues smoothly from 
the phase of dynamical mass generation to the normal phase. Furthermore, the definition (|^) 
of the running coupling is also justified within the conventional framework of Gell-Mann-Low 
renormalization of the gauge model. 

The following two approximations were suggested in ref. f^: 
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Using (|^),and neglecting B{k) in the denominator of (|l|), the following approximate form for 
A{p) in (P may then be derived @] 

A{p) = 1 - f f dk^^{{^)'Q{p -k) + Q{k - p)}. (8) 

Note that the form of the infrared regulator, which has now been introduced explicitly, is- 
foUowing l^-simply a momentum-space cut-off. In [Q an approximate expression for gji was 
given which was based on (H), but with the further approximation of retaining only the term 
with Q{k — p), while neglecting the one involving Q{p — k). This produces the result (always 
forn = 1 in (g)) 
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for < p < e. (9) 



In I, by contrast, the term Q{p — k) in (j^) was retained and the one with Q{k—p) was dropped, 
leading to 

an approximation which was argued to be valid for p ^ a. Whereas g^'^ grows as p decreases, 
g^'^ exhibits the opposite behaviour, and this was interpreted in I as indicating that the true 
dependence of gn on p/a would be significantly slower than that of g^'^ - a situation we 
described in terms of a 'quasi fixed point structure '. 



It is clearly important to establish the extent to which the conclusions drawn from (P) or 
(p!0|) are in fact reliable. Of course, a numerical solution of (|1|) and is always possible (see 
below), but more insight can often be gained from analytical results, provided they are not 
artefacts of a too-crude approximation scheme. In both I and the authors were interested 
in exploring the effects of different ra's in (^). But, as stated above, the case n = 1 seems the 
most relevant. The main purpose of the present note is to exploit the simplifications associated 
with the n = 1 case to investigate more rigorously the issues raised in I and [Q. 



Before embarking on the calculations, we remark that for simplicity we shall present results 
only for one value of N (or ^fo)- We select the number of flavours to lie comfortably in the 
regime where dynamical mass is not generated. From the analysis of and |§] it becomes clear 
that there is a critical number of flavours for dynamical mass generation to occur, and that this 



depends on the infrared cut-off. Although at present we have not studied the corresponding 
critical line in any detail, we can be sure that the value N = 5 will exceed the critical line for the 
regime of infrared cut-offs that are physically interesting to us i.e. 10^^ < 6/a < 10^^. We 
shall present results for the value = 5. 



The simplest improvement of I and 
find (for n = 1) 
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which is simply the obvious combination of (^ and (ID). The quantities g'^' and g'^"' are 
shown in figure 1, from which we see that, as compared with (^, (^) provides essentially 
a smooth approach to the region < p < e where gR is constant (in this approximation). 
In this region, there is a line of non-trivial infrared fixed points, corresponding to the value 
9*R = 5'o/(l + ^ln{e/ a)). By demanding positivity of yl'^i'(p, e) at p = e one can obtain a mild 
restriction on e: 

__3_ 

e > e 90 ct (13) 

which is reasonable for g^ « 1. 



The existence of a non-trivial infrared fixed point was one of the main conjectures of I, relevant 
for the non-fermi liquid behaviour of the system. Actually the same structure is evident in the 
n = 1 curve of figure 1 of 0, though it was not remarked upon there; a similar conjecture was 
also made previously for four-dimensional QCD in ref. |^. The second conjecture made in I, 
also implying a non-fermi liquid behaviour, pertains to a significant slowing down of the rate 
of decrease of g^ with increasing p at intermediate scales of momerd^a, relative to that shown 
by 9r^- Such behaviour is not apparent from figure 1; indeed, g^and g^^ run essentially 
parallel until p is quite close to e. However, the conjectured behaviour at intermediate values 
of p will in fact emerge from the more accurate treatment which now follows. As we shall see, 
the restriction ([l3| ) also disappears, and one gets a non-trivial infrared fixed point structure at 
p = 0. 

In order to improve upon ([TlD , we return to (|l|) with B = 0, n = 1, and a momentum-space 
cut-off e. The running of the coupling gR{p, e) is then defined by (^ with 

^(p,e) = l-^4 r^I{p,k) (14) 
TT'^Jy p-^ Je k 

which we have evaluated numerically. The resulting gR{p,e) is also shown in figure 1, for 
comparison with g^ and gR^- It is indeed apparent that the exactly evaluated gR{p, e) does 
run significantly more sl<2y tha„ either of the approximate quarrtMes J aj;-. thus 
confirming the second conjecture of I. 



In the discussion thus far, we have implemented the infrared regularization by a simple 
momentum-space cut-off e. It may be questioned whether it is consistent, in this case, to consider 
values of the external momentum p lying below e,as we have done in figure 1. Indeed, the 
behaviour of gnip, e) for p near e and below it is, in our view, unlikely to be physical. Detailed 
examination shows that dg^/dp vanishes at p* = e{l + 0[e/a]), where gji reaches a maximum. 
However, for smaller values of the external momenta there is a decrease of g^ towards p = e, so 
that p* is not a fixed point. For values of p less than e, g^ decreases further, leading ultimately 
(it appears) to a fixed point at p = 0. This behaviour seems unphysical, and may signal an 
inconsistency between the vertex ansatz of (^ and the momentum-space cut-off scheme, at small 
momenta. We would like to find an infrared regularization in which both p and the integration 
variable k could range between and a, and which would be expected to lead to smoother 
results, for small p, than those of figure 1. Besides, it is important to check whether the main 
quahtative features of our results are independent of the details of the infrared regularization 
scheme. 



Fortunately, there is a consistent way of implementing the infrared cut-off within the ansatz 
(I), which satisfies these criteria. We keep the limits of integration from to a, and interpret 
the mass function B/A in ([1|) as a (covariant) infrared cut-off in the case of no dynamical mass 
generation. The expression for A then reads 

a \ 1'°' k 

^M-l-^j,ldk^APM, (15) 

and the associated coupling gR^p, 6) is defined by gn^p, S) = go/A{p, 6). This way of introducing 
the infrared (IR) cut-off makes some contact with the finite temperature case |1|, |], where 
the plasmon mass was interpreted as an effective infrared cut-off. The above similarity is, 
however, only indicative. Whether the situation described here carries over intact to the finite- 
temperature regime is at present only an expectation. These are issues that are left open for 
future investigations. 



We may now repeat, for the "5" cut-off case, the calculations decribed above for the "e" 
cut-off. We first consider the use of (|^) in ([l5|) , with the further restriction of keeping only the 
0(fc — p) term, as in 0]. This gives 
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which agrees with (P) for p > e as 5 — 0, and provides a smooth continuation to p = 0. Next, 
we retain both 6 terms, as in (||), leading to 

g^^^'' = go/A(^\p,S) (17) 



where (c.f. (|T|)) 

Again we note that for fixed 6 (|18]) goes smoothly to a constant value as p — > 0. Finally, we 
calculate gR^p, 6) numerically. 



The results are shown in figure 2, which may be compared directly with figure 1. The 
first observation is that the appearance of the three curves is qualitatively very similar in the 
"e" and "5" cases, confirming that the detailed form of the infrared cut-off is not making a 
qualitative difference. Secondly, one sees that g]^^'^ runs parallel to g^^'^ until p is near 6, 
and that both these approximations give significantly faster running than the exact gR{p,S). 
Finally, we stress that there is a non-trivial frixed point at p = 0, which is approached perfectly 
smoothly,in contrast to the unusual behaviour shown in figure 1. The result for gii{p,S) is 
therefore in full agreement with the expectations of I regarding a 'walking technicolour'-like 
behaviour of gn at intermediate scales of momenta. 



To emphasize the insensitivity of these results to the form of the infrared cut-off, we show 
again in figure 3(a) the curves gni^p, e) of figure 1 and gni^p, S) of figure 2, for e = 0.5 and S = 0.5. 
Of course, there is no reason why such a comparison has to be made at the same numerical 
values of e and 6, and we could have adjusted one of them to get closer matching of the results. 
But this would contain no extra physics: figure 3(a) already shows convincingly enough that the 
two crucial elements-the non-trivial fixed point aX p = and the slow running-are independent 
of the form of the cut-off. 



It interesting to note that the above features are characteristic of the presence of a finite 
infrared cut-off. Removal of 6 via 5 — in a smooth manner does not seem to be possible, as 
becomes clear from figure 3(b). From a physical point of view where the infrared cut-off 
is conjectured to be associated with temperature in certain condensed-matter systems whose 
physics the above model is supposed to simulate, this would imply that the above non-trivial 
structure is an exclusive feature of the finite-temperature field theory. 



We add one final comment on the various approximations introduced above. Our exact 
numerical results based on either (|1^) or ( plSj) show that the attractively simple formula (H) (or 
the analogous one with the 6 cut-off), which requires both of the approximations (^ and (|^), 
fails to capture the slow running at intermediate momenta. However, approximation (^ alone 
does appear to be sound, and worth pursuing further. The expression for A based on (|]), with 
the 6 cut-off, is 

^^'nP,5) = l-^- / dk-^^{f{p,k) + f{k,p)} (19) 
n'^N p-^ Jo + 0^ 

where 

2A;3(4aV + 4aV + aV't^ - + a^A;^ + 3apA;4 + 3p2A;4) 
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and we define the corresponding g^ by g^^^ = go/A^^\p,6). The integrals in (|^) can be 
performed analytically, with the result: 
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The interesting point is that gj^ {p, S) turns out to be almost indistinguishable from the exact 
gji{p,6), showing that (P) is actually very reliable. This is shown in figure 4. The availability 
of this very accurate analytical expression for might be useful. For example, taking the 
derivative of (|2l|) with respect to p it is straightforward to show that the non-trivial infrared 
fixed point does lie at p = 0. It is probably of interest to the detailed reader to note that the 
more crude approximation consisting of keeping only the Q{k — p) terms in (|^), which is an 
improved version of what was done in [Q, exhibits a much better agreement with the exact 



answer of fig. 2 and (|21|) than that based on the approximation @. This shows that the 
inadequacy of the approximation of ref. |^ lies mainly in the steps leading to rather than 
in the neglect of modes with momentum below p. 



To conclude, in this short note we have solved the SD equations (P, with the ansatz n = 1 in 
(^), for the wave- function renormalization in the normal phase of the model, where dynamical 
mass is ignored. The solution has been exact, in the sense that it is valid for a wide region of 
momenta < p < a. The solution indicated that the increase of the running coupling is cut-off 
in the infrared, in the form of a non-trivial fixed point, and that there is also a significant 
slowing down of the running of the coupling at intermediate scales, as compared to the case 
of ref. or of Eq. (|11]). This slowing down, or 'walking-technicolour-like' behaviour, has 
been argued ^ to be responsible for the smallness of the coherence length of the dynamical- 
mass-generation phase of the theory. On the other hand, both features, the non-trivial infrared 
structure and the walking-coupling behaviour at intermediate scales, have been argued in I to 
be responsible for deviations from fermi-liquid behaviour, which might have important physical 
consequences, in case the model simulates correctly the physics of the novel high-temperature 
superconductors. We have also argued that the above non-trivial low-energy structure is a 
consequence of a finite infrared cut-off. Thus, from the physical point of view, adopted in I, 
of associating the infrared cut-off with temperature effects, this would imply that the above 



structure is an exclusive feature of the finite-temperature field theory. 
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Figure Captions 



Figure 1: Running coupling constant of large-N resummed QED3 in the normal phase as 
a function of the external momenta p, in the presence of a momentum-space infrared cut-off e. 
In the figure e = 0.1, a = 1, and = 5; the dotted curve shows g^'^ of (y); the top continuous 
curve shows g^^^ of (|Tl|) and ([T^); the bottom continuous curve shows the exact (numerical) 
result, gRip,e). 



Figure 2: As in figure 1, but now employing a covariant infrared-mass cut-off 6 = 0.1. The 
dotted curve is g^^'^ of (|I6|); the top continuous curve is g^^^'^ of (|T^; the bottom continuous 
curve is gRip,5). 



Figure 3: A comparative study of the running of the QED^ coupling in the cases with 
momentum-space (e) and convariant-mass (6) infrared cut-offs for a = 1, A^ = 5. In (a) the 
dotted and continuous curves show gR{p,e = 0.1) and gR{p,S = 0.1) respectively. In (b) is 
shown the value of the coupling gR{p = x,x) at momenta of order of the generic infrared cut- 
off X, as a function of the cut-off (the dotted curve corresponds to the e cut-off, whilst the 
continuous one corresponds to the 5 cut-off). Clearly the cut-off cannot be removed smoothly 
in either case. 



Figure 4: Test of the validity of the analytic result (pi]) , based on the approximation (P), 
in the covariant mass-cut-off 6 case. The running of the coupling based on (|2T|) (dotted curve) 
is compared with the exact numerical result (continuous curve), for a = 1, 6 = 0.1 and A^ = 5. 
The agreement is very good throughout the entire range of external momenta < p < 1. 



